Quantum phases of quadrupolar Fermi gases in optical lattices 
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We introduce a new platform for quantum simulation of many-body systems based on nonspherical 
atoms or molecules with zero dipole moment but possessing a significant value of electric quadrupole 
moment. We consider a quadrupolar Fermi gas trapped in a 2D square optical lattice, and show 
that the peculiar symmetry and broad tunability of the quadrupole-quadrupole interaction results 
in a rich phase diagram encompassing unconventional BCS and charge density wave phases, and 
opens up a perspective to create topological superfluid. Quadrupolar species, such as metastable 
alkaline-earth atoms and homonuclear molecules, are stable against chemical reactions and collapse 
and are readily available in experiment at high densities. 
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Quantum gases of ultracold atoms have provided 
a fresh perspective on strongly-correlated many-body 
states, by establishing a highly tunable environment in 
which both open questions of solid state physics and 
novel, previously unobserved, many-body states can be 
studied [l[ . An important landmark was reached by cool- 
ing and trapping dipolar atoms and molecules, bosonic 
and fermionic [2|, near or into quantum degeneracy, 
which extended the range of features available to quan- 
tum simulation in ultracold atom systems beyond con- 
tact interactions. Numerous exotic states such as su- 
persolids, quantum liquid crystals and bond-order solids 
have been predicted, extended Hubbard models with 3- 
body interactions, and highly tunable lattice spin models 
for quantum magnetism have been proposed 
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crucial feature of the interactions in dipolar gases is their 
anisotropic and long-range character, which is key to the 
intriguing many-body effects that have been predicted. 
In addition, the use of microwave and optical fields to 
shape the symmetry of the interactions between polar 
molecules has been proposed [1, Q . 

In this Letter we propose to study quadrupolar quan- 
tum gases. This constitutes a new class of systems 
in ultracold atom physics, which can be used as a 
platform for quantum simulation of many-body states. 
Quadrupole interactions are most visible for particles 
with zero dipole moment, but which possess a signif- 
icant electric quadrupole moment due to their non- 
spherical shape. The angular dependence of the resulting 
quadrupole-quadrupole interaction is substantially differ- 
ent compared to the dipole-dipole one, due to the higher- 
order symmetry. For atoms or molecules in an optical 
lattice, this allows for a broad tunability of the relative 
strengths and signs of the nearest and next-nearest neigh- 
bor couplings. As a concrete example, below we dis- 
cuss metastable alkaline-earth atoms and homonuclear 



molecules which have comparatively large quadrupolar 
moments. To demonstrate rich many-body effects that 
arise in ensembles of such particles, we derive the quan- 
tum phase diagram of a quadrupolar fermionic gas in an 
optical lattice at half-filling. We find that several uncon- 
ventional phases emerge, such as bond order solids and 
p-wave pairing, and discover the intriguing possibility of 
creating topological ground states of p x + ip y symme- 
try. While dipolar quantum gases were also shown to 
host novel many-body phases, quadrupolar particles are 
available in experiment at higher densities and are stable 
against chemical reactions p| and collapse jjj. 

In order to determine the quadrupole-quadrupole in- 
teraction energy of two atomic or molecular states, we 
consider the potential of a classical quadrupole with mo- 
ment q = J(Sk ■ f — l)dq located at r*o = aligned in 
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FIG. 1. Recipe for realization of quadrupolar particles: 

(a) with alkaline-earth atoms in long-living 3 P2 levels; and 

(b) with homonuclear molecules in rotational states with 
J > 0. (c) Angular "shape" of quadrupolar particles exem- 
plified by |2, M) states. 



2 



^-direction, with k being a unit vector 0, and f = f/r. 
The resulting potential is <fi(f) — ^-(3 cos 2 6—1), where 
cos# = k ■ f. If a second quadrupole with the same 
alignment k is placed at location f, the resulting interac- 
tion energy is V qq = ■ V)(fc ■ V0) - ^(35 cos 4 9 - 
30 cos 2 9+3). The functional form of this potential carries 
over to the quantum description, however, the prefactor 
of this interaction for two | J, M) states has to be obtained 
via a quantum definition of the quadrupole moment q2. 
The latter is a spherical tensor of rank two with com- 
ponents defined as qi, m = - J2k r fcC2,Af(0fc, </> fc ), where 
(rk,9k, 4>k) give the coordinates of the fc-th electron of the 
particle, and C 2 ^i(9 k ,<pk) = Y^W 5 *2,M(0fc> <t>k) are the 
reduced spherical harmonics [8( (here written in atomic 
units). In the angular momentum basis, |J, M), with 
M being the projection of the angular momentum, J, 
onto the quantization axis, the quadrupole operator cou- 
ples the states with A J = 0,±2, so to first order any 
state with J > possesses a nonzero quadrupole mo- 
ment 0. Thus the value of the quadrupole moment can 
be controlled by preparing the particles in a particular 
| J, M)-state, or a combination thereof, using optical or 
microwave fields. In quantum mechanical language, the 
quadrupole interaction is 



(1) 



where (r, 6, 4>) gives the radius- vector between the par- 
ticles, and [q 2 ^ <£> q 2 2 ']4, a is a spherical tensor of rank 
four formed from two quadrupole moment tensors of rank 
two. For both particles prepared in the same | J, M) state, 
Eq. CO) reduces to V qq = V(3 - 30 cos 2 9 + 35 cos 4 9)/r 5 , 
with V = q 2 3{J 2 + J-3Af 2 ) 2 /[4(4J 2 + 4J-3) 2 ], where q 
is defined as q = 2(2, 2|(?2,o|2, 2), which coincides with the 
classical definition [9( . We note that in the classical limit 
of J — > oo, and for M = J, the prefactor V = 3q 2 /16 
of the classical expression is recovered. The quadrupolc- 
quadrupole interaction can be attractive or repulsive de- 
pending on the angle 9. 

Among the particles for which the quadrupolar mo- 
ment is known, the most promising candidates for the ex- 
perimental realization of quadrupolar quantum gases are 
metastable alkaline-earth atoms 0, [T3| and homonuclear 
diatomic molecules Alkaline-earth atoms, such as Sr 
or Yb, can be prepared in metastable st ates, whose 
lifetime exceeds thousands of seconds (il llOj. by optical 



1 In contrast to a dipole which is analogous to a single-headed 
arrow (f) pointing in a particular direction, a quadrupole corre- 
sponds to a double-headed arrow (£) not favoring one direction 
over another, and therefore it can be aligned, but not oriented. 

2 Note the difference with dipole moments which can be nonzero 
only for states of indefinite parity (superposition of odd and even 
J's) 



excitation, cf. Fig.[T](a). Both bosonic and fermionic iso- 
topes of Sr and Yb have been brought to quantum degen- 
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12j. Furthermore, atoms in such metastable states 
can be magnetically trapped (l3| . Ultracold homonu- 
clear molecules, such as Cs2 or Sr2, can be prepared in 
the absolute ground state, 1 S+ (v = 0, J = 0), and then 
transferred to a rotational state with J > 0, using a two- 
photon Raman transition far-detuned from an electron- 
ically excited state [lH , cf. Fig. Q] (b) . While homonu- 
clear molecules are always bosons, fermionic quadrupo- 
lar molecules can be prepared by photo- or magneto- 
association of distinct isotopes of the same species [lij ]. 
For both atoms and molecules the degeneracy of a par- 
ticular J level can be lifted by an external electric or 
magnetic field, F, defining the laboratory z axis, and the 
particles can be prepared in a particular M-component 
possessing a permanent quadrupole moment. The typ- 
ical "shapes" of quadrupolar states are exemplified in 
Fig.Q](c). Both atoms and molecules can be prepared in 
the 1 2, 0) (| 2, 2)) states using two linearly (circularly) po- 
larized photons; the quadrupole-quadrupole interaction 
is equal in this cases and larger than for the |2, 1) states. 

The quadrupole moments for metastable alkaline-earth 
atoms and homonuclear molecules are similarly of the 
order of 10-40 a.u. [ii-fll|. which gives an interaction 
strength, V qq of Eq. ([1]) , on the order of a few Hz at 266 
nm lattice spacing. Furthermore, interactions of the or- 
der of 1 kHz can be achieved for 100 nm lattice spacings 
realizable with atoms trapped in nanoplasmonic struc- 
tures [HI. We note that the dispersion (van der Waals) 
interaction, V dls ~ r~ 6 , is 10 2 — 10 3 times smaller at typi- 
cal optical lattice spacings [l(J , therefore the quadrupole- 
quadrupole interaction dominates the physics of these 
systems. Quantum gases can be trapped for tens of sec- 
onds, so the observation of many-body phases generated 
by these interactions is feasible via the standard tech- 
niques, ranging from time-of-flight detection to noise cor- 
relation and Bragg spectroscopy. 

To illustrate the intriguing many-body effects that can 
arise in these systems, we investigate the quantum phase 
diagram of a system of interacting quadrupolar fermions 
on a square lattice, at half-filling. We assume all parti- 
cles to be prepared in the same state, \J,M), where J 
is the electronic (for atoms) or rotational (for molecules) 
angular momentum, and M is the projection of J on the 
direction F = (9f,4'f) in the laboratory frame given by 
the external field F used to lift the M-degeneracy. The 
particles are confined to a lattice with a lattice constant 
a^, corresponding to the Hamiltonian: 
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t j c\ciCjCj, 
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where t represents the nearest-neighbor hopping and c,; is 
the fermion annihilation operator at the «-th lattice site. 
Throughout the remainder of the paper, we use ql as a 
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FIG. 2. (a) Schematic representation of quadrupolar fermions 
on a square lattice. Alignment of the quadrupoles is given by 
the quantization axis of the external field F, pointing along 
F = (0f,4>f)- The nearest-neighbor interaction is repre- 
sented by green and red solid lines, while the next-nearest 
neighbor interaction is shown in blue, (b) 3D plot showing 
the interactions V& (red), V y (green), and Vx+ y (blue) as a 
function of the angles (9f,4>f); "*" marks the point in the 
vicinity of which both V^,g and V^+y change the sign. 



unit of length and t as a unit of energy. As schemat- 
ically shown in Fig. [5] (a), the interaction strength 
depends on the orientation of the vector connecting the 
quadrupoles, r = r; — r^, relative to the field direction, 
F, via V r = Vij = (ij\V™\ij) = V[3 - 30 (f • F) 2 + 35(f • 
F) 4 ]/r 5 . Thus, one can immediately observe that the in- 
teraction between two quadrupoles can be tuned either 
attractive or repulsive, by changing the orientation of the 
external field F with respect to the interparticle radius- 
vector r. Fig. [2] (b) shows the (Of, 0F)-dependence of 
the interaction matrix element between the nearest- and 
next-nearest neighbors. The richness of the quadrupolar 
interaction becomes apparent in this figure. There are 
several regions in which the signs and the relative mag- 
nitudes of {V x , Vy, V x +y} show distinctive characteristics. 
For example, in the region (9 F < 25°, 0° < <f> F < 45°), 
both nearest- and next-nearest neighbor interactions are 
repulsive, while they all become attractive in the region 
(30° < F < 60°, <f) F ~ 45°). Furthermore, one can iden- 



tify finite regions where either one or two of {Vx, V y , 
V^+g} is attractive while the rest is repulsive. 

Interactions of opposite sign can result in competition 
between quantum phases of different symmetry, resulting 
in frustration. Thus, fermions with dominant quadrupo- 
lar interactions provide an interesting setup for studying 
many-body physics with competing phases. For example, 
in the vicinity of (90°, 45°) both V x and Vg are attractive, 
while V x+y is repulsive (see Fig. [2J . On general grounds, 
one would expect a BCS type ground state resulting from 
condensation of Cooper pairs due to the attractive V x 
and V y couplings. However, the repulsive V±+ y interac- 
tion, if significant, may lead to the insurgence of some 
other phase, and therefore needs to be quantitatively 
accounted for. As another intriguing example, in the 
vicinity of (40°, 5°), V x is strongly attractive while V y is 
strongly repulsive. As we show below, the ground state 
in this region is neither a BCS state nor the conventional 
CDW (such as checkerboard or CDW S ). These two ex- 
amples show that, due to anisotropy, the actual ground 
state may be of an unexpected nature. Exposing the 
true ground state thus demands a theory that is (i) unbi- 
ased with respect to the initial ansatz - condensation of 
particle-particle or particle-hole pairs; and (ii) includes 
fluctuations. 

Issue (ii) can be adequately addressed within the renor- 
malization group (RG) analysis at weak couplings, where 
the low energy physics near the Fermi surface is ex- 
tracted by successively integrating out the high energy 
degrees of freedom [lfj . In order to satisfy criterion (i) 
as well, we employ the exact (or "functional") renor- 
malization group (FRG) which keeps track of all the 
interaction vertices, including both the particle-particle 
and particle-hole channels, and treats all instabilities (in- 



cipient orders) on equal footing [171 ]. Specifically, all 
Feynman diagrams up to one-loop order are included 
for calculating the flow of the most general 4-point ver- 
tex function E/]j(ki, k2, Here k 1: 2,3 are momenta 
on the Fermi surface, with k.4 = ki + k2 — k3, and £ 
the FRG flow parameter starting with the initial value 
I = for the bare quadrupolar interaction. The BCS and 
CDW instabilities are indicated by the diverging flow of 
the corresponding vertex functions C/^(ki, — ki, k2) and 
Ui (ki, — k 2 , ki + Q) respectively Here Q = (±it,±it) 
are the nesting vectors. The most divergent flow indi- 
cates a broken symmetry phase, the order parameter of 
which has an orbital symmetry given by the correspond- 
ing eigenvector of the diverging vertex function. For the 
example considered above, in the vicinity of (90°, 45°), 
the system favors a BCS state with p y orbital symmetry, 
shown in Fig. [3] (b). 

The FRG phase diagram, Fig. [3] (a), features sev- 
eral different kinds of BCS and CDW phases. CDW S 
is the familiar charge density wave phase with a checker- 
board modulation of on-site densities, occurring in re- 
gions where the repulsive interaction between nearest 
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FIG. 3. Quantum phase diagram for quadrupolar fermions on a square lattice, (a) FRG phase diagram in the weak coupling 
limit, V/t = 0.2, and at half filling shown as a function of the magnetic field direction F = (8f,4>f). The point marked with 
"*", where 5 different phases seem to meet, corresponds to V x , Vy, V x + y ~ as shown in Fig. [2](b). (b) The orbital symmetry 
of the CDW and BCS phases shown in (a), obtained from FRG. The p x (p y ) wave CDW has the same orbital structure as the 
Px (p y ) wave BCS. 



neighbors dominates, see Fig. [2 This happens for all 
values of 4>f when Of < 25°, and also for 4>f < 22° 
at large Of > 60°, cf. Fig. [5] In addition, two novel 
types of CDW, CDW Px and CDW Ph , are present. They 
correspond to a checkerboard modulation of the effective 
hopping between nearest neighbors along the x and y di- 
rection respectively, i.e., (cjcj) with r^ — Vj = x or y, with 
the average taken over the many-body ground state. We 
refer to these phases as to bond order solids (BOS). In 
comparison, the s-wave CDW order corresponds to mod- 
ulations of (cfci). We refer to these phases as to bond 
order solids (BOS). Furthermore, we find a small region 
of CDW s +d that involves a mixture of extended s- and 
d- waves. The individual contributions of the extended 
s-wave and d-waves vary smoothly within the CDW s +<j 
region. Together they give rise to a checkerboard modula- 
tion of effective hopping between the next-nearest neigh- 
bor sites. The CDW P;E , CBW Py , and CDW s+d can be 
thought of as a 2D generalization of the bond-order-wave 
phase occurring in the extended Hubbard model in one 
dimension [HI]. BOS was recently shown to arise in 2D 
for dipolar fermions 0]. Compared to dipolar systems, 
BOS occupies a significantly larger region of the param- 
eter space for quadrupolar interactions (e.g., it is sta- 
bilized as soon as 9p approaches 25°). Moreover, the 
angular dependence of quadrupolar interactions is sub- 
stantially more complex, resulting in two BOS phases of 
p y and p x symmetries, appearing at small and large Of, 
respectively. Interestingly, these two phases occur in the 
regions where V x and V y are comparable in magnitude 
but opposite in sign, i.e., CDW Pa .(p ) is stabilized when 
Vx( y ) is repulsive while Vg( x ) is attractive, cf. Fig. [2](b). 



Thus, quadrupolar Fermi gases are well suited for ex- 
ploring the properties of nonzero angular momentum (i.e. 
unconventional) density wave phases. 

The remainder of the phase diagram is occupied by 
two BCS phases. They roughly occur in regions where 
both Vx and V y are attractive. Our FRG analysis shows 
that the BCS phase can be stable even through the next- 
nearest neighbor interaction is weakly repulsive. We find 
that the symmetry of the BCS order parameter is p x or 
p v , depending on whether V x or V y is more attractive. 
Along the line of Op ~ 65° these two BCS phases are 
degenerate. This raises the possibility of realizing p x +ip y 
topological superfluid order using quad rupolar gases. By 
analogy with the proposal of Ref. [19( , using an ac field 
to periodically modulate the direction of (0f,4>f), one 
can lift the degeneracy and engineer the chiral p x + ip y 
state. 

In conclusion, we have shown that ultracold Fermi 
gases with quadrupole-quadrupole interactions can be 
used to study unconventional BCS, CDW, and topologi- 
cal phases, and gain insight into the physics of compet- 
ing ground states. While we have focused on the specific 
case of a square lattice at half-filling, the functional RG 
methods of this work can be applied to study other fill- 
ings and lattice geometries. Temperatures achieved for 
degenerate Fermi gases of alkaline-earth atoms in experi- 
ment are T = 0.26 T F and T = 0.37 T F respectively 0]. 
The optimal T c for the CDW and BCS phases predicted 
here is estimated to be on the order of 0.03 Tp, for inter- 
mediate couplings, V ~ t. Thus these many-body phases 
seem to be within experimental reach in the near future. 

We note that quadrupole-quadrupole interactions are 
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on average less anisotropic compared to the dipole-dipolc 
ones, which results in an ultracold gas more stable against 
collapse and paves the way to study many-body physics 
in three spatial dimensions. Furthermore, quadrupolar 
interactions occur in numerous subfields of physics, from 
molecular photofragmentation [2(| and structure of /- 
electron compounds 21 to nuclear reactions [22} and 
gravitation of black holes [13]. Thus, potential appli- 
cations of the proposed quantum simulation platform 
spreads far beyond the many-body physics of fermionic 
gases. Finally, we note that ground-state atoms can be 
provided with significant quadrupolc moments by means 
of Rydberg dressing, i.e. admixing a highly-excited elec- 
tronic state possessing a large quadrupolc moment with 
far-detuned laser light [24 1. 
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